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We investigate the dynamical structure of advective accretion flow around 
stationary as well as rotating black holes. For a suitable choice of input pa- 



in 
in 

. rameters, such as, accretion rate (M) and angular momentum (A), global 

accretion solution may include a shock wave. The post shock flow is located 



at few tens of Schwarzchild radius and it is generally very hot and dense. This 
successfully mimics the so called Compton cloud which is believed to be re- 
sponsible for emitting hard radiations. Due to the radiative loss, a signiflcant 
energy from the accreting matter is removed and the shock moves forward 
towards the black hole in order to maintain the pressure balance across it. 
We identify the effective area of the parameter space {Ai — A) which allows 
accretion flows to have some energy dissipation at the shock (A£). As the 
dissipation is increased, the parameter space is reduced and flnally disap- 
pears when the dissipation is reached its critical value. The dissipation has a 
profound effect on the dynamics of post-shock flow. By moving forward, an 
unstable shock whose oscillation causes Quasi-Periodic Oscillations (QPOs) 
in the emitted radiation, will produce oscillations of high frequency. Such an 
evolution of QPOs has been observed in several black hole candidates during 
their outbursts. 
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1 INTRODUCTION 

In a significant work on the pro spect of sliock formation in an accretion disk around a black 
hole, Chakrabarti &: DaJ ( 2004 ) showed that in order to have a stable shock, the viscous dis- 
sipation inside a flow must have an upper limit, beyon d whi ch the Ran kine-Hugoniot condi- 



tions cannot be satisfied. In 



Das fc Chakrabarti! (120041 ) and iDad (120071 ) , the bremsstrahlung 



and synchrotron cooling were also added to dissipate away the heat generated from viscous 
dissipation. However, it is well known that the post-sh ock region emits the hard X-rays in 



a black hole candidate ( Chakrabarti fc Titarchuk 



19951 ) and some amount of energy is lost 



through radiation from the post-shock region. This radiative loss primarily comes from the 
thermal energy of the flow and takes place via thermal Comptonization. In a self-consistent 
shock condition, this radiative loss must also be incorporated. In the present paper, we 
quantitatively show how the energy loss at the shock affects the location of the shock itself 
around stationary as well as rotating black holes. As energy dissipation is increased, the 
post-shock flow pressure gets reduced causing the shock front to come closer to the black 
hole in order to maintain the pressure balance across it. Accordingly, the dynamical prop- 
erties of standing shock waves would directly be related to the amount of energy discharge 



ition, the mass out 


low rate which is believed to be gener- 


( Chakrabarti 


1999; 


Das et al. 


2001; 


Das & Chattopadhvav 



20081 ). would also be affected by the energy discharge at the shock location. Therefore, it 
is pertinent to understand the response of energy dissipation on the formation of standing 
shock wave. In this paper, we precisely do this. 

The plan of our paper is the following: in the next Section, we present the equations 
governing the flow and the procedure we adopted to solve these equations. In §3, we show how 
the Rankine-Hugoniot conditions at the shocks must be modified when energy dissipation is 
present. In §4, we show the results of our computations. Finally in §5, we present concluding 
remarks. 
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2 GOVERNING EQUATIONS AND SONIC POINT ANALYSIS 

We start with a steady, thin, rotating, axisymmetric, accretion flow around black holes. We 
assume smaller accretion rates, so that the flow radiatively ineflicient and behaves essentially 
as a constant energy flow as in Chakrabarti (1989). We assume a polytropic equation of state 
for the accreting matter, P = Kp'^ , where, P and p are the isotropic pressure and the matter 
density, respectively, 7 is the adiabatic exponent considered to be constant throughout the 
flow and is a constant which measures the entropy of the flow and can change only at 
the shock. Since we ignore viscous dissipation the angular momentum of the flow A = x^q 
is also constant everywhere. However, we assume that the main dissipation is concentrated 
in the immediate vicinity of the post-shock flow, which would be the case if the thermal 
Comptonization is the dominant process. The flow height is determined from the condition of 
being in equilibrium in a direction perpendicular to the equatorial plane. Flow equations are 
made dimensionless considering unit of length, time and the mass as GMbh/c^, GMbh/c^ 
and Mbh respectively, where G is the gravitational constant, Mbh is the mass of the black 
hole and c is the velocity of light. 

In the steady state, the dimensionless energy equation at the disk equatorial plane is 
given by (jChakrabartilll989l ). 



1 .0 



2 7 



(1) 



where, S is the speciflc energy, is the radial velocity and a is the adiabatic sound speed 
defined as a - 



I^P/ p. Here, effective potential due to black hole is denoted by $. In the 
present study, the effect of gravity is taken care of by two different potentials. To repre- 
sent Schwarzschild black hole, we use Paczyhski-Wiita (jPaczynski fc Wiital 119801 ) potential 
i^ pw) and for the Kerr b l ack h ole, we consider pseudo-Kerr potential {^pk) introduced 



by 



Chakrabarti fc Mondall (120061 ). It has been adequately shown that these potentials accu- 



rately mimic not only the geometry of the space-time, but also the dynamics of the flow. In 
fact the error for n ot using full general relativist ic treatment has been shown to be at the 
most a few percent ( jChakrabarti fc Mondalll2006h . The expressions for Paczyhski-Wiita and 
pseudo-Kerr effective potential are respectively given by. 



1 



PW 



2x2 2{R - 1) 



and 
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PK 



B + - AAC 
2A 



where, 



A 



2x2 ' 



B = -1 + — + 



X 



C= 1 



1 



R — Xq X ' 

Here, x and R represent the cyhndrical and spherical radial distance from the black hole 
when the black hole itself is considered to be located at the origin of the coordinate system. 
Here, Xq = (0.04+0.97afc+0.085a^)/2, oo = 2ak/{x^+alx+2al) and = (^x2_2^^(j^2^y(^^2^ 
a\ + 2a\/ x), a is the red shift factor, represents the black hole rotation parameter defined 
as the specific spin angular momentum of the black hole. 

The mass flux conservation equation in the steady state apart from the geometric factor 
is given by, 

M = dpxh{x), (2) 
where, M is the mass ac cretion rate consid ered to be constant, and h{x) represents the 



half-thickness of the flow (IChakrabarti 



19891 ) which is expressed as. 



h{x) 



X 



7$ 



(3) 



R 



Here, $ 



R 



\dR 



z«x 



and z is the vertical height in the cylindrical co-ordinate system where 



R = a/x2 + 2;2. By using the polytropic equation of state and the definition of the adiabatic 



sound speed, we get the equation for the so-called 'entropy accretion rate' (IChakrabarti 
19891) as: 



M = 



7$ 



(4) 



R 



where, v = (7 + l)/(7 — 1). 

In order to form shocks, accretion flow must be supersonic at some point, i.e., stationary 
flow must pass through a son ic point. We deriv e the sonic point conditions following the 



standard sonic point analysis (IChakrabarti 



19891 ). We, therefore, differentiate Eqs. (1) and 



(4) with respect to x and eliminate terms that involve derivatives of a to obtain radial 



Studies of dissipative standing shock waves around black holes 5 
velocity gradient which is expressed as, 



(M_N 
dx~D' 

where, the numerator N and the denominator D are respectively given by, 

d^' 



(5) 



x(7 + l) dx (7 + 1)$^ dx ^ ' 



and 



D^d- , (56) 
(7 + 1)^? 

Here, a subscript 'e' signifies that the quantity was calculated on the equatorial plane. 

In order to have the flow smooth everywhere, the numerator and the denominator in 
Eq. (5) must vanish simultaneously. These considerations allow us to obtain the sonic point 
conditions. Therefore, setting D = 0, we have. 



^l^^c) = j:^^al{x,) M,= ^-^. (6a) 

The subscript 'c' signifies that the quantity was computed at the sonic points. Note that 
the Mach number [M(xc)] at the sonic point is not unity as in the other models since the 
vertical equilibrium model is considered. The other condition comes while setting N — 
and is given by, 

-1 

(66) 



«c(^c) = (7 + l)''^'^' 



1 d^ji 



X dx 



Sound speed at the sonic point can now easily be calculated from Eq. (6b) and it must be 
always positive. Flow may possess at most three sonic points outside the black hole horizon. 
Among them, the closest one from the black hole is called inner sonic point whereas the 
furthest one is known as outer sonic point. The nature of the sonic point solely depends on 
the value of velocity gradient (Eq. 5) at the sonic point. Usually, d'd/dx has two values at 
the sonic point: one is for accretion fiow and the other is for winds. A sonic point is referred 
as saddle-type when both the derivatives are real and of opposite sign. When derivatives are 
real and of same sign, the sonic point is called as nodal type. If the derivatives are complex, 
the sonic point is center-type. The saddle type sonic point receives special importance as 
transonic flow can only pass through it. Moreover, for standing shock transition, flow must 
have multiple saddle type sonic points. We solve Eqs.(l-6) to obtain a full set of global flow 
solutions where standing shock may be present. 
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3 SHOCK CONDITIONS AND SHOCK INVARIANT 

In order to have a shock, the flow velocity must jump discontinuously from super-sonic to 
sub-sonic branch. This discontinuous shock jump is characterized by four flow variables, 
namely, shock location (x^), radial velocity (-(9), sound speed (a) and entropy constant {K) 
respectively and are given by. 



X = {7 a) 

A^ = Mxs)-^-{xs)i (76) 

Aa = a+{xs) — a_(xs), (7c) 

and 

AK = K+- K^, {Id) 

where, the subscripts '-' and denote the quantities before and after the shock. Here, we are 
dealing with a quasi two dimensional flow where pressure and density are averaged over in the 



vertical direction. In this case, the Rankine-Hugoniot shock conditions (ILandau fc Lifshitz 



19591 ) are given by. 



8+ = 8-- A8, (8a) 
(II) the self-consistent pressure condition, 

W+ + = W. + (86) 

and (III) the baryon number conservation equation, 

M+ = M_. (8c) 



Here , W and S denote the vertically averaged pressure and density of the flow (IMatsumoto et al 



mi)- 



Because of our choice of radiatively inefficient flow, i.e., flow with low accretion rates, 
the energies 8^ and 8^ are nearly constant. However, in presence of dissipation at the shock, 
A8 will be non-zero. Since the dissipation is expected to be mostl y through thermal Comp 



toniz ation, which may cool the flow from ~ lOOkeV to ~ 3keV fiChakrabarti fc Titarchuk 



19951 ). we assume that the energy dissipation in the post-shock region reduces the tempera- 



ture of the flow and the loss of energy is proportional to the temperature difference between 
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the post-shock and the pre-shock flows similar to the well known Kirchoff's law, i.e., 

AS = AS'{al-al), 

where, AS' is the proportionality constant which will be used as a parameter. In reality, this 
is a function of the electron an d low energy photon contents o : 



be calculated in any situation (jChakrabarti fc TitarchukI Il995l ). Since the thermal energy 



the accretion flow and can 



content (the final term in Eq. 1) can be anywhere between ten to thirty percent of the rest 
mass energy, a few percent of thermal energy can be easily extracted. Below, we show that 
maximum dissipation for a black hole of a ~ 0.8 is about 6.5%. 
(I) the energy conservation equation. 



Using the 



at the shock (IChakrabarti 



governmg ec 



natio ns and shock conditions, we compute an invariant relation 



19901 ) which is given by. 



C 



[M+(37-l) + 



2 12 



[M_(37-l) + 



2 12 



2 + (7 - l)Ml + 2(7 - 1)AS' 2 + (7 - l)Ml + 2(7 - 1)AS' 
where, M's are the Mach numbers of the flow. 

Simplifying Eq. (9) we have the following relations. 



2X 



and 



where, 



A' = (37-i)'-c(7-i), y 



(37-l)'-C(7-l)' 



(9) 



(10) 



(11) 



and Z = 4. 



2(37 - 1) - C - C(7 - 1)A^ 

The flow will pass through the shock if the shock invariant condition is satisfied (Eq. 9) 
at some point between the two saddle type sonic points. 



4 RESULTS 

Fig. 1 shows a phase space trajectory of accretion flow around a Schwarzschild black hole 
where the logarithmic radial distance is plotted along the horizontal axis and the Mach 
number (M = i^/a) is plotted along the vertical axis. Arrows indicate the direction of the 
accreting flow. Here, the flow parameters are A = 3.56, Ai = 0.08 x 10^^ and 7 = 4/3. 
Sub-sonic accreting flow passes through the outer sonic point (O) and becomes super-sonic. 
Depending on the energy dissipation at the shock (AS), shock conditions are satisfied at some 
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log(x) 

Figure 1. Phase space diagram of accretion flow in the M-log(x) plane. Flow pajrameters are A = 3.56 and M. = 0.08 x 10~^. 
The energy at outer sonic point is obtained as £ = 0.13452 x 10~^. The vertical lines indicate the possible shock transitions 
Xs = 27.75 (big-dashed-dot) , Xs = 17.02 (small-dashed-dot) and Xs = 10.001 (dotted) corresponding to the energy dissipation 
at the shock Af = 0.699 x lO'^, A£ = 0.37 x lO'^ and AS = 0.985 x 10-^ respectively Shocks form closer to the black hole 
horizon as the energy dissipation at the shock is enhanced. 




log(xJ 

Figure 2. Variation of the energy dissipation at the shock as a function of shock location, (a) Accretion rate of the flow is 
chosen as M = 0.01 x 10~®. Curves are plotted for a set of angular momentum starting form A = 3.46 (left most curve) to 
A = 3.82 (right most curve) in the interval of d\ = 0.04. (b) Accretion rate of the flow is chosen as M = 0.01 x 10~®. The 
leftmost curve is drawn for A = 3.14 and other curves are for increasing angular momentum with an interval of dX = 0.04. 
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Figure 3. Plot of the energy dissipation at the shock with the shock location for a given angular momentum, (a) A = 3.70 
and (b) A = 3.22. See text for more details. 

particular location and the flow makes a discontinuous jump (shock-transition) from super- 
sonic branch to sub-sonic branch there. After the shock transition, the flow momentarily 
slows down and subsequently picks up it velocity to cross the inner sonic point (1) before 
entering into the black hole. Since the flow losses a part of its energy at the shock in the form 
of thermal energy, the pressure in the post-shock region is diminished. As a result, the shock 
moves forward to maintain a pressure balance across it. The possible shock transitions at Xs 
are shown as the vertical lines for various energy dissipations at the shock. For a given set 
of input parameters (A = 3.56, = 0.08 x 10~^), the shock forms closer to the black hole 
horizon as the energy dissipation at the shock increases. For instance, Xs = 27.75, AS = 
0.699 X lO"-'^ (big-dashed-dotted), = 17.02, AS = 0.37 x lO'^ (small-dashed-dotted) and 
Xs = 10.0, AS = 0.985 X 10"^ (dotted) respectively 

In Fig. 2(a-b), we present the variation of shock location in the logarithmic scale as a 
function of the energy dissipation at the shock. In the upper panel, we show the results 
obtained by using PW potential while in the lower panel, we use the PK potential. In 
Fig. 2a, the accretion rate is chosen as = 0.01 x 10^^ and adiabatic index is 7 = 4/3. 
Different curves are drawn for a set of angular momentum starting from A = 3.46 (left most 
curve) to A = 3.82 (right most curve) in the interval of dX — 0.04. In Fig. 2b, we consider 
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TVI = 0.01 X 10^^ and black hole rotation parameter = 0.5. In this plot, the angular 
momentum is varied from A = 3.14 (left most curve) to A = 3.30 (right most curve) in the 
interval of d\ = 0.04. Both the Figures indicate that the shock forms closer to the black 
hole horizon as the energy dissipation increases monotonically for a given input parameters 
(A^,A). The conclusion is valid for both the rotating and non-rotating black holes. Note 
that there is a cut-off in each individual curve. At very high energy dissipation RHCs are 
not satisfied and no steady shocks can form. However, non-steady shocks may still form as 
accreting matter continues to possess multipl e sonic points w here entropy at the inner sonic 



point is higher than the outer sonic point (jPas et al.l 1200 ll ). In this case, shock location 



becomes imaginary and shock starts oscillating with a time period which is comparable to 
the inverse of observed quasi periodic oscillation (QPO) frequencies of black hole candidates. 
In order to investigate this physical process explicitly, a rigorous time dependent analysis 
is required which is beyond the scope of our present work. As the angular momentum is 
increased, maximum energy dissipation at the shock first increases, becomes maximum at 
some A and then decreases. For a non-rotating black hole, the maximum energy dissipation 
at the shock (ASmax) is around 3.5% for accretion flows with angular momentum A = 3.70 
and for rotating black hole AS^ax is around 4.7% for A = 3.22. 

A diagram similar to Fig. 2(a-b) is drawn in Fig. 3(a-b) for various entropy accretion 
rates keeping the specific angular momentum fixed. In the upper panel, we present the 
results for flows with angular momentum A = 3.70 accreting onto a non- rotating black 
hole. Here, M is varied from M = 0.01 x 10"^ (leftmost) to 0.21 x 10^^ (rightmost) in 
the interval of dAi = 0.01 x 10~^. In the lower panel, we present a similar plot for flows 
of angular momentum A = 3.22 accreting around a rotating black hole with the rotation 
parameter = 0.5. The accretion rate is chosen to be 7V1 = 0.01 x 10~^ (left most curve) 
to 0.16 X 10^^ (right most curve) in the interval of dAi = 0.05 x 10^®. In both the Figures, 
for a given accretion rate, the shock front moves forward with the enhancement of energy 
dissipation at the shock. Interestingly, there is a significant reduction in maximum energy 
dissipation (ASmax) at the shock as Ai is increased. Thus the parameter space for standing 
shock is reduced with the increase of the energy dissipation at the shock. 

In Fig. 4, we classify the parameter space spanned by the accretion rate (Ai) and the 
angular momentum (A) on the region which allows standing shocks for various values of 
energy dissipation at the shock. In the left and the right panels, the parameter spaces are 
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Figure 4. Modification of the region of the parameter space which forms a standing shock as a function of energy dissipation 
(A£) at the shock. Left paneh For non-rotating black hole. Right panel: For rotating black hole (a^ = 0.5). Parameter space 
shrinks with the increase of A£. 



obtained for non-rotating and rotating black holes. Here, rotating black hole has ak = 0.5. 
In the left Figure, the parameter space under the solid curve is obtained for a dissipation- 
free accretion flow. The region under the dotted and the dashed curves are obtained for 
AS = 1.5 X 10~^ and 3.0 x 10~^ respectively. Similarly, in the right panel, the parameter 
spaces separated by solid, dotted and dashed curves are obtained for AS — 0.0, 1.0 x 10~^ 
and 4.0 x 10~^. These figures show that the effective region of the parameter space shrinks 
with the increase of the energy dissipation at the shock in the lower angular momentum 
side. Above a critical limit of AS, this region disappears completely. 

An important part of understanding an accretion flow around a rotating black hole is 
to compute the maximum energy dissipation (AS^ax) at the shock as a function of the 
black hole rotation parameter a^. In Fig. 5, we present the variation of AS^ax with a^. For 
instance, for an weakly rotating black hole, AS^ax is around 2.5%. For increasing a^, AS^ax 
increases and reaches upto 6.5% for a rapidly rotating black hole. However, the maximum 
available energy is around 5.7% and 13% of the rest mass energy for = and = 0.8 
respectively. Thus we find that maximum released energy at the shock is about 40 — 50% of 
the available energy. 
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Figure 5. Variation of A£max with black hole rotation parameter Ofc. See text for more details. 



5 CONCLUSIONS 



It is well known that the accretion flows around a stellar mass bl ack hole emit hard X-rays 
from a Compton cloud located close to the inner edge of the disk (jChakrabarti fc Titarchuk 
19951 ). This Compton cloud is the post-shock region which also produces jets and outflows 
from the disk. The Comptonization process removes mainly the thermal energy of the flow. In 
this paper, we show that at standing shock the maximum release of energy could be as high 
as 40 — 50% of the maximum available energy. Since the radiative loss removes a signiflcant 
thermal pressure from the inflow, the location of the shock is affected. In particular, the 
shock seems to form closer to the black hole when the dissipation is increased due to the 
drop of the post-shock pressure which no longer can support the shock in its original place. 
The implication of this could be profound. First of all, the observed QPOs are considered 
to be due to the oscillation of this shock. During outbursts, the QPOs of several black 
hole candidates, such as, XTE J1550-564 and GRO J1655-40 evolve rapidly. The frequency 
changes from several mHz to a few tens of Hz in the rising phase of the outbursts. Since the 
QPO frequency is inversely proportional to the infall time of the post-shock flow, such an 
evolution indicates that the shocks actually propagated towards the black hole at the rising 
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phase of the outbursts (IChakrabarti et al.ll2008l . 120091 ) . Our mechanism discussed here gives 
a plausible explanation of such a phenomenon. 
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